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OBJECTIVE — To quantitatively assess the strength and shape of the association between 
blood 25-hydroxy vitamin D [25(OH)D] levels and incident risk of type 2 diabetes. 

RESEARCH DESIGN AND METHODS— A systematic search of the MEDLINE and 
Embase databases and a hand search of references from original reports were conducted up 
to 31 October 2012. Prospective observational studies that assessed the association between 
blood levels of 25(OH)D and risk of incident type 2 diabetes were included for meta-analysis. 
DerSimonian and Laird's random-effects model was used. A quadratic spline regression analysis 
was used to examine the shape of the association with a generalized least-squares trend test 
performed for the dose-response relation. 

RESULTS — A total of 21 prospective studies involving 76,220 participants and 4,996 incident 
type 2 diabetes cases were included for meta-analysis. Comparing the highest to the lowest 
category of 25(OH)D levels, the summary relative risk for type 2 diabetes was 0.62 (95% CI 
0.54-0.70). A spline regression model showed that higher 25(OH)D levels were monotonically 
associated with a lower diabetes risk. This inverse association did not differ by sex, duration of 
follow-up, study sample size, diabetes diagnostic criteria, or 25(OH)D assay method. A linear 
trend analysis showed that each 10 nmol/L increment in 25(OH)D levels was associated with a 
4% lower risk of type 2 diabetes (95% CI 3-6; P for linear trend < 0.0001). 

CONCLUSIONS — Our meta-analysis showed an inverse and significant association between 
circulating 25(OH)D levels and risk of type 2 diabetes across a broad range of blood 25(OH)D 
levels in diverse populations. 



Low vitamin D status is prevalent in 
many populations and has become a 
common public health problem 
worldwide (1,2). Vitamin D is well-known 
for its essential role in calcium homeostasis 
and bone health (2). Emerging evidence 
from both in vitro and in vivo studies has 
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suggested extraskeletal effects of vitamin D, 
including on insulin action and secretion 
(2,3). Population studies have provided 
further support to the hypothesis that low 
vitamin D status, as assessed by circulat- 
ing 25-hydroxy vitamin D [25(OH)D] 
levels, is associated with impaired (3-cell 



function, insulin resistance, and impaired 
glucose intolerance and thereby may be 
associated with higher risk of type 2 di- 
abetes (4-7). An inverse association be- 
tween 25(OH)D levels and prevalent type 
2 diabetes has been shown (4,8-10); 
however, a temporal relationship cannot 
be established from such cross-sectional 
studies that are subject to bias due to the 
possibility of reverse causation. Recently, 
several prospective observational studies 
have reported a significant association be- 
tween high circulating levels of 25(OH)D 
and lower incidence of type 2 diabetes (11- 
17). However, no association was observed 
in other studies (18-21). In addition, indi- 
vidual studies have been underpowered to 
examine this relation across a broad range 
of circulating 25(OH)D levels. As high- 
lighted by the recent report on vitamin D 
from the Institute of Medicine (IOM) 
(1,22), available evidence on the relation 
between vitamin D status and type 2 diabe- 
tes remains inconsistent and inconclusive, 
and there is a need for further research to 
clarify the optimal levels of 25(OH)D for 
nonskeletal outcomes, including type 2 
diabetes, and to assess whether there is 
a nonlinear relationship between 25(OH)D 
and diabetes risk. Assessment of possible 
threshold levels of 25(OH)D will greatly 
advance our understanding of the magni- 
tude and shape of the association of vita- 
min D with incidence of type 2 diabetes. 

Available evidence from prospective 
studies with circulating 25(OH)D levels 
measured in baseline blood samples col- 
lected before the onset of disease allow us to 
address the temporal relation of vitamin D 
status with risk of type 2 diabetes. We 
therefore conducted a meta-analysis of 
prospective studies from various popula- 
tions to quantify the association between 
circulating 25(OH)D levels and subsequent 
risk of type 2 diabetes and examine possible 
threshold effect of 25(OH)D levels. 

RESEARCH DESIGN AND 
METHODS 

Data source and searches 

Relevant studies were identified by search- 
ing MEDLINE and Embase databases for 
all published articles up to 31 October 
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2012, using the search terms "vitamin D," 
"25-hydroxy vitamin D," "25(OH)D," "1,25- 
dihydroxy vitamin D," "l,25(OH)2D," "cal- 
cidiol," "calcitriol," "insulin resistance," 
"insulin secretion," "insulin sensitivity," "glu- 
cose intolerance," "glucose metabolism," 
"beta-cell function," "type 2 diabetes," and 
"diabetes." The search was further re- 
stricted to English-language articles, hu- 
man studies, and adult subjects aged > 19 
years. Additional studies were retrieved 
through a hand search of references from 
original reports. All prospective studies on 
this topic were considered eligible if they 
provided data on the relationship between 
baseline circulating levels of 25(OH)D and 
risk of type 2 diabetes. 

In the first round of screening (n = 
542), 470 articles were excluded for at 
least one of the following reasons: studies 
that did not study circulating 25(OH)D as 
an exposure or type 2 diabetes as an out- 
come (200 articles), studies among chil- 
dren and adolescence populations (4 
articles), nonhuman studies (15 articles, 
including chemistry, animal, cell line, and 
isolated tissue studies), reviews/editori- 
als/guidelines/letters/commentaries (185 
articles), cross-sectional or retrospective 
studies (61 articles), and case reports (5 
articles) (Fig. 1). In the second round of 
screening, full-text articles were retrieved 
(n = 72). Sixty-three articles were further 
excluded for at least one of the following 
reasons: updated data available from 
other studies in the same population or 
duplicate publications (3 articles), 



diabetes complications as major out- 
comes (33 articles), association measures 
for diabetes-related risk factors (19 articles), 
and lack of directly measured 25(OH)D lev- 
els (2 articles using predicted score). Our 
search strategy and inclusion/exclusion cri- 
teria resulted in a total of 2 1 independent 
prospective studies (extracting from 15 ar- 
ticles) being included in the current meta- 
analysis. 

Data extraction 

Three investigators (Y.S., L.W., and L.C.D.G.) 
independently reviewed each eligible 
article and extracted relevant data exam- 
ining the prospective associations of cir- 
culating levels of 25(OH)D with type 2 
diabetes risk. Differences, if any, were 
reconciled through group discussion. 
Data extracted include population source, 
study design, follow-up period, sample 
size, subject characteristics (age, sex, and 
comorbid conditions), 25(OH)D assay 
methods, diabetes end points, and main 
study findings. When results were avail- 
able only on different subpopulations in 
the same study (13,14,23), we separately 
extracted data for each subpopulation 
(study characteristics described in Supple- 
mentary Table 1). 

Statistical analysis 

To provide quantitative evidence from all 
studies and maximize statistical power for 
hypothesis testing, we performed meta- 
analyses using DerSimonian and Laird's 
random-effects model with inverse-variance 



(SE) weighting of individual study results 
(24) when data could be combined. The 
summary relative risks (RR) with 95% CI 
were calculated for the highest versus the 
lowest category of 25(OH)D levels. We 
also performed subgroup meta-analyses 
to explore potential effect modification 
by prespecified factors, including sex, 
follow-up duration, sample size of cases, 
adjustment for BMI, adjustment for BMI 
and other metabolic parameters (includ- 
ing hypertension, lipids, or inflammatory 
biomarkers), and 25(OH)D assay meth- 
ods. Heterogeneity across studies was as- 
sessed by Cochran Q statistic. Because 
the Cochran Q statistic has low statistical 
power (25), we also calculated the I 2 and 
H statistics to reflect between-study het- 
erogeneity. The percentage of I 2 —25 
(J 2 = 25), 50 (J 2 = 50), and 75 (I 2 = 75) 
indicates low, medium, and high hetero- 
geneity, respectively (26). The H statistic 
is a complement to the I 2 statistic in assess- 
ing study heterogeneity; an H <1.2 indi- 
cates little heterogeneity, and an H >1.5 
raises caution regarding notable hetero- 
geneity (26). We visually assessed publi- 
cation bias using Begg modified funnel 
plots, in which the RR was plotted on a 
logarithmic scale (log[RR]) against its SE 
from each study. Publication bias was 
also assessed by two formal tests: Begg 
adjusted rank correlation test and Egger 
regression asymmetry test (27,28). The 
former was used to examine if there was 
significant correlation between the effect 
estimates and their variances, and the latter 
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Figure 1 — Flow chart of study selection. 
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uses an inverse-variance weighted regres- 
sion of the effect sizes on their precision 
(the inverse of SE) to test whether the in- 
tercept deviates significantly from zero. 

We tested for possible nonlinear as- 
sociation between circulating 25(OH)D 
levels and type 2 diabetes using a quadratic 
polynomial spline regression analysis, 
considering a possibility of data perturba- 
tion at an extreme end of the exposure 
range (29). We also used the method de- 
scribed by Greenland and Longnecker (30) 
for the dose-response analysis to compute 
the linear trend from the correlated RRs 
and 95% CIs across categories of 25(OH)D. 
The median level of 25(OH)D in each cat- 
egory was assigned to the corresponding 
RR when reported in the study. If not re- 
ported, the values assigned were the means 
or the midpoint of the lower and upper 
bound in each category. For extreme 
open-ended categories, half the width of 
the adjacent category was subtracted (for 
the lowest category) or added (for the up- 
permost category) to obtain the midpoint. 
The numbers of case and control subjects 
or person-years by category were also col- 
lected if available. 

All analyses were performed using the 
STATA statistical software (Version 10.1; 
STATA Corp., College Station, TX). Sta- 
tistical significance was defined as two- 
tailed a < 0.05. 

RESULTS— We included 15 publica- 
tions (11-16,18,19,23,31-33) that met 
our inclusion criteria in our meta-analysis, 
which provided data from 2 1 independent 
studies [11 cohort studies (11-13,17,18, 
20,21,31-33), 8 nested case-control stud- 
ies (14,15,19,23), and 2 case-cohort stud- 
ies (16,33)] (Fig. 1). These prospective 
studies comprised a total of 76,220 partic- 
ipants and 4,996 incident diabetes cases. As 
summarized in Supplementary Table 1, 13 
studies were population-based (12- 
14,16,17,20,21,23,32,33), and 2 studies 
included postmenopausal women only 
(18,19). Nineteen studies included largely 
whites (11-18,20,21,23,32,33), and 2 
studies included multiple racial/ethnic 
populations (19,31). The outcome of dia- 
betes was ascertained using a combination 
of criteria, including diabetes-specific 
pharmacotherapy, diabetes-related hos- 
pitalization, self-report, and glycemic 
status information (11-16,18,19). Five 
studies used the American Diabetes Asso- 
ciation glycemic criteria for diagnosis 
(20,21,23,31,33). The assay method for 
25(OH)D varied across studies. Radio- 
immunoassay was used in 8 studies 



(14,15,18,21,33), chemiluminescent im- 
munoassay was used in 10 studies (11- 
13,16,17,19,23,32), and only 3 studies 
used liquid chromatography-tandem 
mass spectrometry, which is the gold stan- 
dard to measure 25(OH)D (20,31,33). 

Hypothesis testing 

Fig. 2 shows the summary RR of type 2 
diabetes comparing the highest to lowest 
category of 25(OH)D levels for hypothe- 
sis testing. The summary RR was 0.62 
(95% CI 0.54-0.70), indicating a signifi- 
cant inverse association between baseline 
25(OH)D levels and risk of type 2 diabe- 
tes. The P value for the Cochran Q test 
was 0.21, the H statistic was 1.1 (1.0- 
1.5), and the I 2 statistic was 19% (0- 
53), indicating no evidence of significant 
between-study heterogeneity. The Begg 
funnel plot for the visual assessment of 
publication bias showed that smaller 
RRs with large SEs tended to be above 
the horizontal line, indicating a possibil- 
ity of publication bias in favor of small 
studies with null findings (data not 
shown). However, the Egger test (P = 
0.39) and the Begg test (P = 0.99) did 
not show evidence of publication bias. 

Sensitivity analysis 

We conducted a sensitivity analysis to 
assess the extent to which individual 
studies with extremely large RRs influ- 
enced the summary RR. The exclusion of 
the study by Anderson et al. (11) that 
included the largest sample size and RR 
estimate reduced between-study hetero- 
geneity but did not appreciably change 
the summary RR (0.64; 95% CI 0.56- 
0.73; P for Q statistic = 0.27, H = 1.1 
[1.0-1.4], and I 2 = 15% [0-50]). Omit- 
ting four studies from the earliest publi- 
cation by Knekt et al. (14), which 
included small sample size and heteroge- 
neous results, also did not substantially 
influence the summary RR (0.61; 95% 
CI 0.55-0.68; P for Q statistic = 0.46, 
H = 1.0 [1.0-1.4], and I 2 = 0% [0-51]). 
The most recent study with the smallest 
number of cases (n = 37) yielded large 
variance in the effect estimate; removing 
this study led to almost the same sum- 
mary RR (0.62; 95% CI 0.56-0.70; P 
for Q statistic = 0.35, H = 1.0 [1.0-1.3], and 
I 2 = 9% [0-44]). 

Sources of heterogeneity 

We evaluated potential effect modifiers in 
stratified analyses (Table 1). The inverse 
association between 25(OH)D and diabe- 
tes risk was consistently observed in men, 



women, and a mixed population of men 
and women. Associations tended to be 
stronger for men than for women alone, 
but the difference did not reach statistical 
significance. Nor did duration of follow- 
up, study sample size, diabetes diagnostic 
criteria, and 25(OH)D assay methods in- 
fluence the summary RRs. Several studies 
presented RRs for 25(OH)D levels and 
type 2 diabetes with and without adjust- 
ment for metabolic variables, including 
BMI, hypertension, markers of glycemia 
and insulin sensitivity, plasma lipid lev- 
els, or inflammatory markers. The associ- 
ations between 25(OH)D levels and 
diabetes risk were attenuated but re- 
mained statistically significant after ad- 
justment for BMI and hypertension and/ 
or other biomarkers: the RRs were 0.52 
(95% CI 0.46-0.58) and 0.63 (0.56- 
0.72) in the models without and with ad- 
justment for these covariates, respectively 
(P for interaction = 0.02). Likewise, ad- 
justment for BMI alone slightly attenuated 
the association (P for interaction = 0.06). 

Dose-response analysis 

Fig. 3 shows the dose-response relation 
between circulating 25(OH)D and risk of 
type 2 diabetes. Overall, the test for a lin- 
ear relation across the range of 25(OH)D 
from 20 up to 160 nmol/L was signifi- 
cant (P < 0.0001). The RR for type 2 di- 
abetes was 0.96 (95% CI 0.94-0.97) per 
10 nmol/L increment in 25(OH)D. A sig- 
nificantly lower risk of type 2 diabetes 
became evident when 25(OH)D approxi- 
mated 50 nmol/L. Because most studies 
had 25(OH)D ranges < 100 nmol/L, cur- 
rent evidence for the relation of 25(OH)D 
of > 100 nmol/L with type 2 diabetes was 
weak. 

CONCLUSIONS— On the basis of 21 
prospective studies including 4,996 in- 
cident cases of type 2 diabetes and 76,220 
nondiabetic controls, our meta-analysis 
showed a significantly inverse association 
between 25(OH)D levels and incidence of 
type 2 diabetes. The association was 
consistent and did not differ appreciably 
by sex, study size, follow-up duration, 
diabetes diagnosis criteria, or 25(OH)D 
assay methods. The association was at- 
tenuated but remained significant after 
adjustment for BMI and/or intermediate 
biomarkers. 

Blood 25(OH)D levels, which reflect 
all sources of vitamin D exposure and 
have a half-life of 2 to 3 weeks, have been 
widely used as a surrogate of vitamin D 
status (2,3). Due to differences in assay 
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Figure 2 — A random-effects meta-analysis of 21 independent prospective studies with adjusted RR and 95% CI of type 2 diabetes in relation to serum 
25(OH)D levels (the highest category versus the lowest category). AusDiab, Australian Diabetes, Obesity and Lifestyle Study; DPP, Diabetes 
Prevention Program; Ely, Medical Research Council Ely Study; EPIC, European Prospective Investigation into Cancer-Norfolk Study; EMC, Finnish 
Mobile Clinic Health Examination Survey; Hoorn, Hoorn Study; Inter99, Inter99 Study; MFH, Mini-Finland Health Survey; MONICA1, Danish 
MONICA1 survey ; MONICA/KORA, Monitoring of Trends and Determinants in Cardiovascular Disease/Cooperative Health Research in the Region 
of Augsburg Study; NHS, Nurses' Health Study; PS, Pizarra Study; SDPP, Stockholm Diabetes Prevention Program. 



methods and population characteristics, 
there is no consensus on the cutoff values 
defining vitamin D insufficiency or defi- 
ciency. Vitamin D insufficiency has been 
previously reported to range from levels 
of 40-75 nmol/L (16-30 ng/mL), and vi- 
tamin D deficiency is generally defined as 
levels of <50 nmol/L (20 ng/mL) (2,34). 
Using 25(OH)D levels <50 nmol/L 
(20 ng/mL) to define vitamin D deficiency, 
data from 4,495 adults in the National 
Health and Nutrition Examination Survey 
2005-2006 showed that the overall prev- 
alence of vitamin D deficiency was 41.6%, 
with the highest rate found in blacks 
(82.1%), followed by Hispanics (62.9%), 
others (57.6%), and whites (30.9%) (35). 
The IOM's new guidelines for vitamin D 
intake and desirable blood levels, which 
were based on a systematic scientific review 
of available evidence at the time (1,22), 
state that a level of 50 nmol/L (20 ng/mL) 
of serum 25(OH)D is needed to maintain 
bone health for most individuals. The IOM 



report called for additional research to clar- 
ify the relation of vitamin D levels with 
nonskeletal outcomes (1,22). 

The optimal 25(OH)D levels for type 
2 diabetes prevention remain unknown. 
Epidemiological evidence relating lower 
circulating vitamin D levels to hypergly- 
cemia, insulin resistance, or type 2 diabetes 
primarily derives from cross-sectional re- 
ports (4,8-10). Prospective data are lim- 
ited and have been inconclusive (36). Our 
dose-response curve showed an inverse 
and significant relation between 2 5 (OH) 
D and type 2 diabetes across a broad 
range of 25(OH)D levels in diverse popu- 
lations. Our results also confirmed that 
baseline 25(OH)D levels >50 nmol/L 
were significantly associated with a lower 
risk of type 2 diabetes, although further 
studies with higher power are required 
to provide more stable estimates of this 
association. 

The observed inverse association 
did not differ by sex, study sample size, 



duration of follow-up, diabetes diagnosis 
criteria, and 25(OH)D assay methods. An 
early pooled analysis of 2 nested case- 
control studies of 412 incident cases of 
type 2 diabetes and 986 control subjects 
found a strong inverse association in men 
but not in women (14). However, this find- 
ing of sex-specific relation of 25(OH)D 
with type 2 diabetes was not confirmed 
by subsequently published studies nor 
by our meta-analysis. It is likely that 
such sex difference resulted from residual 
confounding or statistical fluctuation due 
to small sample size. Adiposity is another 
important covariate in vitamin D-diabetes 
relation. The interrelationship between vi- 
tamin D status and adiposity is complex 
and possibly bidirectional. Increased stor- 
age of 25(OH)D in adipose tissue and less 
sun exposure due to limited mobility and/ 
or excessive subcutaneous fat deposits in 
obese individuals will lead to low circulat- 
ing levels of 25(OH)D (2,3). Conversely, 
vitamin D may directly affect adiposity 



care, diabetesj ournals . org 



Diabetes Care, volume 36, May 2013 



1425 



Vitamin D levels and type 2 diabetes 

Table 1 — Subgroup analyses for the relation between circulating 25(OH)D and type 2 diabetes 
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26 (0-65) 




Diagnostic tools for type 2 diabetes 












0.49 


OGTT 


7 


0.65 (0.53-0.80) 


0.41 


1.0 (1.0-1.9) 


2 (0-71) 




Self-report 


4 


0.71 (0.52-0.96) 


0.20 


1.2 (1.0-2.1) 


35 (0-77) 




Others 


10 


0.58 (0.48-0.70) 


0.24 


1.1 (1.0-1.6) 


22 (0-62) 





OGTT, oral glucose tolerance test. 



and other metabolic parameters such as 
dyslipidemia, hypertension, and systemic 
inflammation that mediate the pathway 
from vitamin D status to type 2 diabetes 



(36,37). In this sense, adiposity could be 
considered as a confounder or an inter- 
mediate variable. Adjustment for adipos- 
ity and other obesity-related metabolic 



parameters may be an overadjustment, 
possibly underestimating the true associa- 
tion. Finally, with the exception of one 
study (31), all studies used a single mea- 
surement of 25(OH)D at baseline, which 
is not a time-integrated measure of vita- 
min D status. 25(OH)D was also assayed 
using different methods in previous stud- 
ies. There was substantial intra- and inter- 
assay variation in 25(OH)D levels (38). 
However, the summary RR of type 2 di- 
abetes in relation to 25(OH)D did not dif- 
fer by assay methods. 

Our findings of an inverse relation 
between 25(OH)D and type 2 diabetes are 
supported by biological evidence that 
vitamin D may be implicated in the 
pathogenesis of diabetes and its compli- 
cations (36,37). A large body of literature 
has suggested that optimal vitamin D ho- 
meostasis is essential for insulin action 
and secretion (2,3), two fundamental fea- 
tures in the pathogenesis of type 2 diabe- 
tes. Clearly, direct evidence from ongoing 
(39) and future clinical trials of higher- 
dose vitamin D supplementation is war- 
ranted to clarify any beneficial effects of 
vitamin D on primary prevention of type 
2 diabetes. 
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Figure 3 — Relation between the risk of type 2 diabetes and baseline levels of 25(OH)D in 18 in- 
dependent prospective studies included in the meta-analysis. The relation is modeled by the quadratic 
spline regression. Circles indicate RR in each study. The circle size is proportional to the precision 
of the RR (inverse of variance). The gray -shaded region shows the 95% CIs around the regression 
line. 
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Findings from our meta-analysis of 
observational studies will augment and 
complement findings from randomized 
trials of the effect of vitamin D supple- 
ments on type 2 diabetes. Although ran- 
domized trials are critical for establishing 
cause-and-effect relationships between 
vitamin D supplementation and health 
outcomes, they will not address all the 
potential questions because of their fixed 
dose (or at most, a few doses) of vitamin D 
and narrow range of 25(OH)D levels 
achieved by supplementation. Prospec- 
tive observational studies allow us to 
assess 25(OH)D thresholds for diabetes 
risk across a broad spectrum of 25(OH)D 
levels. The largest randomized trial of 
vitamin D supplement to date, the Wom- 
en's Health Initiative (WHI) Clinical Trial, 
has evaluated vitamin D plus calcium 
supplementation for fracture prevention 
in >36, 000 postmenopausal women 
(40). Secondary analysis of the WHI trial 
with 33,951 initially nondiabetic women 
did not observe any effect from daily in- 
take of 1,000 mg elemental calcium plus 
400 IU vitamin D3 on incident diabetes 
over 7 years of follow-up (40). Of note, a 
dose of 400 IU/day vitamin D in the WHI 
raised median levels of serum 25(OH)D 
from 42.3 to only 54.1 nmol/L (-12 
nmol/L), which is below the optimal value 
of 75 nmol/L or more for skeletal and 
nonskeletal health including type 2 dia- 
betes (41). Similarly, secondary analysis 
of another large randomized trial, the 
Randomized Evaluation of Calcium Or 
vitamin D (RECORD) trial, did not observe 
any effect from daily intake of 800 IU vi- 
tamin D3 on incident diabetes over 2-5 
years, although such a daily dose raised 
serum 25(OH)D from 38 to 62 nmol/L 
on average (42). 

Our meta-analysis has some limita- 
tions. First, the observational nature of 
prospective studies included in our anal- 
ysis cannot rule out residual confound- 
ing, although the consistency of our 
results across multiple strata and sensi- 
tivity analyses minimizes the likelihood 
that residual confounding explains the 
findings. Second, as in any meta-analysis, 
publication bias is possible, although our 
visual examination and formal tests did 
not suggest the presence of substantial 
publication bias. Finally, limited data 
from existing prospective studies lacked 
sufficient power to detect potential sex- or 
ethnicity-specific 25(OH)D thresholds or 
dose-response relations. 

In conclusion, our meta-analysis of 
21 prospective studies demonstrates a 



significant inverse association between 
circulating 25(OH)D levels and risk of 
incident type 2 diabetes in a dose- 
response manner across a wide spectrum 
of 25(OH)D levels. Direct evidence from 
future randomized trials is warranted to 
clarify a cause-and-effect relationship be- 
tween vitamin D and type 2 diabetes. 
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